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Abstract 

We investigate a two-dimensional classical A-vector model with a nonlinear inter- 
action (1 + <jj • (Tj) p in the large- A limit. As observed for N = 3 by Blote et al. 
[Phys. Rev. Lett. 88, 047203 (2002)], we find a first-order transition for p > p c and 
no finite-temperature phase transitions for p < p c . For p > p c , both phases have short- 
range order, the correlation length showing a finite discontinuity at the transition. For 
p = p c , there is a peculiar transition, where the spin-spin correlation length is finite 
while the energy-energy correlation length diverges. 



PACS: 75.10.Hk, 05.50.+q, 64.60.Cn, 64.60.Fr 



The two-dimensional Heisenberg model has been the object of extensive studies which 
mainly focused on the 0(iV)-symmetric Hamiltonian 



H 



(ij) 



' J- 



where cr^ is an iV-dimensional unit spin and the sum is extended over all lattice nearest 
neighbors. The behavior of this system in two dimensions is well understood. It is disordered 
for all finite (3 [0] and it is described for (3 — > oo by the perturbative renormalization group 
The square-lattice model has been extensively studied numerically 0-0, checking 
the perturbative predictions [PHI5| and the nonperturbative constants |16|-[18|. 

In this paper we study a more general Hamiltonian on the square lattice; more precisely, 
we consider 

H = -N/3 W(l + cr x ■ tr x+(l ), (2) 

where W(x) is a generic function such that W(2) > W(x) for all < x < 2, in order to 
guarantee that the system orders ferromagnetically for (3 — > oo. A particular case of the 
Hamiltonian (0) has been extensively studied in the years, the case in which W(x) is a 
second-order polynomial. Such a choice of W(x) gives rise to the so-called mixed O(N)- 
K.P model |19|-281, which is relevant for liquid crystals [29-34J and for some orientational 
transitions |5H . 

rP + b and found 



In a recent Letter [36 



the authors analyzed a model with W(x) = ax p 
an additional first-order transition for p large enough. Here, we will study the same model, 
finding an analogous result: for p > p c « 4.537857 a first-order transition appears, the 
correlation length — and in general, all thermodynamic quantities — showing a finite discon- 
tinuity. Note that the appearance of a first-order transition in nonlinear models is not a 

oo it was already shown in Ref. 



new phenomenon. Indeed, for TV 
transition appears in mixed 0(A^)-MP Ar 
of interest to understand the behavior for p 



20 



that a first-order 
models for some values of the couplings. It is 





found a 



p c . For such value of p, Ref. 
peculiar phase transition: while the spin-spin correlation length remains finite, the energy- 
energy correlation length diverges. Here, we will show that the same phenomenon occurs 
for iV = oo. However, at variance with what observed in Ref. ||36|1 , the critical theory shows 
mean-field — not Ising — behavior. 

Let us consider the Hamiltonian (||D on a hypercubic d-dimensional lattice. We normalize 
W(x) by requiring W'(2) = 1 so that in the spin-wave limit 



H 



N/3 



dx d^cr ■ d^cr. 



(3) 



We also fix W(l) = so that H = for a random configuration. Then, we introduce two 
new fields X Xfl and p Xfl in order to linearize the dependence of the Hamiltonian on the spin 



coupling. We write 



exv[N(3W(l+a x -(T x+fl )] 



dp X fj,d\xfj, 



cxp 



N(3 



\ xfl (1 + <J X ■ <T 



X + fl 



-p x ») + N(3W{p 



Xfl ) 



(4) 



2 



As usual in the large-iV expansion, we also introduce a field fi x in order to eliminate the 

N{3 



constraint a 2 , = 1. Thus, we write 



-to fal ~ 1) 



5 (crl - l) ~ J dfx x exp 
With these transformations we can rewrite the partition function as 

Z = J Y\[dPxixd\ XIJ ] Y\l d tod(7x} e 



(5) 



MA 



(6) 



where 



XfJ, 



We perform a saddle-point integration by writing 



A 



X/.1 

Pxfi 

to 



a + X XfJ ,, 
r + Pxf,, 
1 + to- 



A standard calculation gives the following saddle-point equations 



d/3(l - r) + - \(2d + m 2 )I(ml) - ll = 0, 
a - 2W'(t) = 0, 







- -Km 2 



2 a 

2\ 



(9) 



where we set 7 = a(2<i + mg)/2 



(27r) d p 2 + ml 



2 • 



(10) 



and p 2 = 4 sin 2 p M /2. The variable m has a simple interpretation: it is related to the 
spin-spin correlation length by £ CT = l/m . From Eq. (§) we obtain finally 



where 



Jlln 



t = r(m ) = 2 + — £ 



2ci 2dl(ml)' 
The corresponding free energy can be written as 

F = -pdW(r) + ~ log I {ml) + ^L(m 2 ) 



(12) 



(13) 
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Figure 1: Function (3 (mo) = I(mo)/W'(r) vs mo, for p = 4, 4.5, 5, and 5.5. For any p, 
P(mo) — ► oo for mo — > 0. 



^K)= |^log(p 2 + ^). (14) 



where 

L(m ° } = / (2*)' 

Focusing now on the two-dimensional case, let us show that, for any W(x), the spin-spin 
correlation length is always finite, i.e. £ CT = oo, i.e. m = 0, only for /3 — oo. Note first that 
r = 2 (resp. t = 1) for m = (resp. m = oo) and that r(m ) is a strictly decreasing 
function of mo- Thus, W'(r) is finite for all mo- Then, since 1(0) = +oo, we find that 
£ CT = oo only if /3 = oo, i.e. £ CT is finite for all finite (3. 

We want now to discuss the behavior for (3 — > oo. From Eq. (]TTj) , we see that /3 — > oo 
for m — > and possibly for m — > mj, where W'(r(mi)) = 0. If there is more than one 
solution, the relevant one corresponds to the lowest free energy. Now, for (3 — * oo, we can 
simply write |38| F m -2f3W(r). Since r(0) = 2 and W(2) > W(t) for all < r < 2 
because of the ferromagnetic condition, the relevant solution is the one with mo — > 0. Then, 
using 

I (mo) = -— log— J + O(m 2 \ogml) (15) 



2tt ° 32 



for m — > 0, we obtain 



= 3 2e - 2 ^+^"( 2 )/ 2 [l + O^" 1 )], (16) 
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Figure 2: The free energy F(f3) for p — 5. There is a critical point C for f3 c ~ 1.543. 



in agreement with the standard perturbative renormalization-group predictions [39 . 

Let us now discuss the possibility of first-order phase transitions, which may arise from 
the presence of multiple solutions to Eq. (O). As in Ref. P5| , we consider 



W(x) 



2 /x\p 
p \2, 



2 i-p 
p 



(17) 



In Fig. |T] we report the function (3 (mo) = J(mo)/W'(r), for p = 4,4.5,5,5.5. For p = 4,4.5, 
for each (3 there is a unique solution m and thus there are no phase transitions. On the 
other hand, for p = 5, 5.5 there is the possibility of multiple solutions, in which case the most 
relevant is the one that gives the lowest free energy. For p = 5, we report the free energy in 
Fig. |[ We observe a first-order transition for (3 ~ 1.543 with a finite discontinuity of the 
correlation length, A£ CT w 16.2, and of all thermodynamic quantities. A numerical analysis 
of the gap equation (|TT|) shows that a first-order transition exists for all p > p c ~ 4.537857. 
For p = p c , the thermodynamic functions are nonanalytic for (3 = (3 C ~ 1.33472. In this case, 

- p e « -0.035726(m - m 0c ) 3 + O[(m - m 0c ) 4 ], (18) 



where m 0c ~ 0.387537. Consequently, repeating the discussion of Ref. 

£ a {p) w 2.5804 + 7.8682(/3-/5 c ) 1/3 + ---, 



(19) 
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E((3) « 0.162274 + 0.314385(/3- ) 5 c ) 1/3 + (20) 
C(fi) « 0.104795(/5-/5 c )" 2/3 + ---, (21) 

where and C are respectively the energy and the specific heat per site. Note that C(/3) 
diverges at the critical point, indicating that, although spin-spin correlations are not critical, 
criticality is observed for energy-energy correlations. Indeed, consider 

D Q (k) = e lHx - y) (Q(l + ct x ■ <r x+tl ); Q(l + a y ■ <r y+u )), (22) 

where Q(x) is an arbitrary regular function. For N —>■ oc, 

D Q (k) = [Q , (T)] 2 J2(M-k);Mk)), (23) 

so that 

ND ^ = {w^) c ^- (24) 

It follows Dq(0) ~ (P — Pc)~ 2 ^ 3 for any function Q(x). Thus, all correlation functions of the 
energy show a critical behavior. In order to compute the associated correlation length, we 
determine Dq(k) for arbitrary k. We obtain 

Nn (k) = 2[Q f (T)]*[A 2 (k)Mk)-A 1 (k)*\ 

Q[ ' P^Wir^Aoi^-pWir^A^A^-A^kyy 1 ] 

where 



a ns f d 2 Q (E M cosg„ 

A n (k) = / — ^ J ^ . (26) 

J W ^fc/ 2 ) + m%) [{q^kj2) + mgj 

For (3 — > fi c and k — > 0, we have 

DQ(k)- 1 = a(J3 - p c ) 2/3 + bk 2 + 0(k 4 ) } (27) 

with a, b 7^ 0. Thus, the energy-energy correlation length £e(P) behaves as 

~ (/3 - /3 c r 1/3 , (28) 



i.e. z/g = 1/3. We thus confirm the results of Ref. |36[ on the existence of the critical theory 
for p = p c , although we disagree on the nature of the critical behavior. Indeed, Ref. p6| 



suggested a = 1 — 1/5, with 5 assuming the Ising value 5 = 15. Instead, we find the mean- 
field value 5 = 3. It is unclear how our large- N result is compatible with what observed for 
N = 3. Indeed, the universality argument of Ref. |$(| would predict Ising behavior for any 
value of N. This issue deserves further investigations. 

We thank Henk Blote and Henk Hilhorst for many useful comments. 
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